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Abstract 

Recent theoretical developments in using the Interacting Boson Model to 
describe nuclear structure effects in fusion reactions below the Coulomb barrier 
are reviewed. Methods dealing with linear and all orders coupling between the 
nuclear excitations and the translational motion are discussed, and the latter is 
found to lead to a better description of the barrier distribution data. A system¬ 
atic study of the available data (cross sections, barrier and spin distributions) 
in rare-earth nuclei is presented. 


1 Introduction 

It is by now well established that to describe fusion cross sections one needs to include 
coupling of relative motion to other degrees of freedom 10. These extra degrees 
of freedom yield a distribution of barriers || and consequently enhance the cross 
section below the barrier ||, |5| . Enhancement of fusion cross sections due to coupling 
of levels in colliding nuclei to the relative motion, has opened a new avenue for testing 
nuclear structure models. The recent experimental determination of average angular 
momenta 0 and the barrier distributions 0 have provided even more stringent tests 
for such models. 

The natural language to study multidimensional barrier penetration is the cou¬ 
pled channels formalism, which was widely used in investigating subbarrier fusion 
phenomena [l|. An alternative formulation of the multidimensional quantum tun¬ 
neling is given by the path integral formalism [|J. An algebraic nuclear structure 
model significantly simplifies evaluation of the path integral. The Interacting Boson 
Model (IBM) of Arirna and Iachello [§] is one such model which has been successfully 
employed to describe the properties of low-lying collective states in medium heavy 


1 


nuclei. Here attempts to use IBM in describing nuclear structure effects in subbarrier 
fusion are reviewed. 

Path integral formulation of this problem, as sketched in the next section, requires 
analytic solutions for the nuclear wave functions. In our first attempt in using IBM 
to describe nuclear structure effects in subbarrier fusion, we employed the SU(3) limit 
of IBM 0 . However, the SU(3) limit corresponds to a rigid nucleus with a particular 
quadrupole deformation and no hexadecapole deformation, a situation which is not 
realized in most deformed nuclei. Thus analytic solutions away from the limiting 
symmetries of the IBM are needed for realistic calculations of subbarrier fusion cross 
sections. In a parallel development, a 1/N expansion was investigated [|T0| for the IBM 
which provided analytic solutions for a general Hamiltonian with arbitrary kinds of 
bosons. This technique proved useful in a variety of nuclear structure problems where 
direct numerical calculations are prohibitively difficult. As we briefly describe in the 
next section, using the 1/N expansion technique in the path integral formulation 
of the fusion problem [[id]] makes it possible to go away from the three symmetry 
limits of IBM, in particular, arbitrary quadrupole and hexadecapole couplings can be 
introduced. 

Recently a series of recent high-precision measurements were carried out at the 
Australian National University (ANU) where distributions of fusion barriers were 
determined directly from the fusion data IBM based analysis of fusion cross 
sections describes the ANU data well ii |T3| , |T4| , especially when the higher 
order coupling effects are included. 

In the next section, first the influence functional method is summarized and its 
application to the linear coupling of translational motion to the structure of target 
nuclei is described. The significance of nonlinear couplings is discussed next, and a 
Green’s function method is introduced to describe such nonlinear effects. Section 3 
includes attempts to describe subbarrier fusion data with these techniques. 


2 Algebraic Models in Subbarrier Fusion 

The Hamiltonian for the multi-dimensional quantum tunneling problem relevant to 
subbarrier fusion is 

H = ~V 2 + H 0 (() + r,£) (1) 

where r is the relative coordinate of the target and projectile, and £ represents any 
internal degrees of freedom of the target. // 0 (£) and H in t (r, £) describe the internal 
structure of the target nucleus and its coupling to the relative motion, respectively. 
In Eq. (1), the propagator to go from an initial state characterized by the relative 
radial coordinate r* (the magnitude of r) and internal quantum numbers n* to a final 
state characterized by 77 and rif may be written as 

K(r f: nf,T;r i: ni, 0 ) = J V[r(t)}e lS ^’ T)/h (nf\U in t(r(t),T)\rii), 
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where S(r,T) is the action for the translational motion and U int satisfies the differ¬ 
ential equation 


[H 0 + H int }U int , (3) 

1. (4) 

We want to consider the case where r\ and ry are on opposite sides of the barrier. In 
the limit when the initial and final states are far away from the barrier, the transition 
amplitude is given by the S'-matrix element, which can be expressed in terms of the 
propagator as fl|] 
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where pi and pf are the classical momenta associated with 7y and rj. In heavy ion 
fusion we are interested in the transition probability in which the internal system 
emerges in any final state. For the tth partial wave, this is 

u=D 5 »/.«.( £ )r> (e) 

n f 

which becomes, upon substituting Eqs. (2) and (5), 
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Here we have assumed that the energy dissipated to the internal system is small 
compared to the total energy and taken pf outside the sum over final states. We have 
also identified the two-time influence functional in Eq. (7) as 

p M {r(t),f;r(t),T) = (ni\U? nt (f(t), T)U int (r(t), T)\m), (8) 

where the completeness of final states was used. Eq. (8) displays the utility of the 
influence functional method when the internal system has symmetry properties. If 
the Hamiltonian in Eq. (3) has a dynamical or spectrum generating symmetry, i.e. 
if it can be written in terms of the Casimir operators or generators of a given Lie 
algebra, then the solution of Eq. (3) is an element of the corresponding Lie group [|J. 
Consequently the two time influence functional of Eq. (8) is simply a diagonal group 
matrix element for the lowest-weight state and it can be evaluated using standard 
group-theoretical methods. 





If the Interacting Boson Model is used to describe the effects of nuclear structure 
on fusion, H 0 = // IB m, and there are two choices for H int . One possibility is to take 
H int to be of the form of the most general one-body transition operator for IBM 

tfint = V(r) + E«*ii(0[6}6i] W • H^(r), (9) 

kjl 


where bj and 6/ denote the boson creation and annihilation operators (bo = s, b 2 = d, 
etc). The k sum runs over k — 2,4,.. . 2Z max , with Z max representing the maximum spin 
of boson operators. Odd values of k are excluded as a consequence of the reflection 
symmetry of the nuclear shape and the k = 0 term is already included in the bare 
potential V(r). The form factors otkji(r) represent the spatial dependence of the 
coupling between the intrinsic and translational motions. The interaction term given 
in Eq. (9) is an element of the SU(6) algebra for the original form of the IBM with 
s and d bosons, and is an element of the SU(T5) algebra when q bosons are included 
as well 

To simplify the calculation of the influence functional, we can perform a rotation 
at each instant to a frame in which the z-axis points along the direction of relative 
motion. Neglecting the resulting centrifugal and Coriolis terms in this rotating frame 
is equivalent to ignoring the angular dependence of the original Hamiltonian. In this 
approximation, the coupling form factors become independent of i and only m = 0 
magnetic sub-states of the target are excited [|Hj. For heavy systems, the neglected 
centrifugal and Coriolis forces are small. We take the scattering to be in the x-y 
plane. Then making a rotation through the Euler angles b = (0, 7r/2, 0), we can 
write the Hamiltonian as the rotation of a simpler Hamiltonian depending only on 
the magnitude of r 

H = R(b)H^(r)R\b). (10) 

Since in Eq. (1) H 0 and H *. are rotationally invariant, H mt is the only term whose 
form is affected by this transformation. Hence we introduce the rotated interaction 
Hamiltonian H^(r), given by 
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4 > jlm('^') 
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If we assume now that the form factors otkji(f ) are all proportional to the same 
function of r, then the Hamiltonian commutes with itself at different times and 
hence we can write the two-time influence functional as 

/ | i dtH^ 0 } (r(t))/h —i f T dtH^°hr(t))/h\ \ f-\ A\ 

Pm = \ni\e J ° mtV e J o m tv uw (14) 


Since the exponents of the two operators in the influence functional commute, Pm 
becomes the matrix element of an SU(6) transformation between SU(6) basis states, 






in other words it is a representation matrix element for this group and can be easily 
calculated using standard techniques. The two-time influence functional for the sd- 
version of IBM was calculated in Ref. [[□]] and, for the particular case of the SU(3) 

limit, in Ref. ||. 

The above method corresponds to a linear coupling approximation. Due to the 
exponential nature of fusion cross sections, it is not clear at the outset whether one can 
include the effects of the higher order terms in H- mt through a renormalization of the 
coupling strength. We next discuss an alternative method, whereby one can include 
the effects of coupling to all orders in the interaction Hamiltonian. In this case the 
solution of Eq. (3) is no longer an element of the appropriate group. Nevertheless, 
we can exploit the symmetry properties of the resolvent operator directly without 
utilizing its path integral representation. Such a Green’s function approach has also 
been used to study quantum tunneling in a heat bath 
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To include the effects of couplings to all orders, the interaction Hamiltonian in 
Eq. (1) is written as 

#int(r,0 = Vcoui(r,0 + Kuc(r,0- (15) 

Here the Coulomb part is 


Vcoui(r,0 = 


Z\Z 2 e 2 


Z x Z 2 e 2 


D3»: 


(r > Ri), 
(r < Ri), 


(16) 


and the nuclear part is taken to have the Woods-Saxon form, 


Kuc(r, 0 = -V 0 1 + exp 


r-Ro-RMr^) 


-l 


(17) 


In Eqs. (16) and (17), Ro is the sum of the target and projectile radii and R\ is the 
mean radius of the deformed target. O is a general coupling operator between the 
internal coordinates and the relative motion which we take to be 




(18) 


The Vk represent the strengths of the various multipole transitions in the target nu¬ 
cleus. In the standard IBM with s and d bosons, the possible transition operators 
have k = 2,4; k = 0 and odd-fc values being excluded for reasons stated after Eq. (9). 
The quadrupole and hexadecapole operators are given by 


T (2) = [s^d, + + x[d^d\^ 2 \ 

T ( 4) = [d f d] (4) . 


(19) 

( 20 ) 


We adopt the “consistent-Q” formalism of Casten and Warner |17j, in which x hi 
Eq. (19) is taken to be the same as in R IB m (fitted to reproduce the energy level 
scheme and the electromagnetic transition rates of the target nucleus), and thus is 
not a free parameter. 















To calculate the fusion cross section to all orders, we consider the resolvent oper¬ 
ator for the system 


G + (E) = ---(21) 

The basic idea is to find the interaction representation which diagonalizes H int . Note 
that, in general, it is not possible to diagonalize Hi BM and H int simultaneously, so the 
effect of the excitation energies in the target nucleus are ignored in this approach. 
To this end, we need to identify the unitary transformation which diagonalizes the 
operator O 

O d = UOU\ (22) 

so that we can calculate its eigenvalues and eigenfunctions 


O d \n) = h n \n). 


(23) 


Assuming the completeness of these eigenfunctions 

H|n)("l = C 


(24) 


one can write the matrix element of the resolvent as 
{Z f ,r f \G + {E)\£ i ,r i ) = ^ f ,r f \U^[E + -H k -H int (r,d)]~ 1 U^ i ,r i ), 

= (0; r 7l wt [ E+ ~ H k - H int (r,6 d )\ 1 ^2\n)(n\U\^ i ,r i ), 

n 

= £«/|WV>(»M«<r/|G:|r ( >, (25) 


where 


Gi(E) = 


E* - H t 


(26) 


Hint (/’, h n ) 

G+(E) in Eq. (26) is the resolvent operator for one-dimensional motion in the poten¬ 
tial H- mt (r,h n ), whose fusion cross section can be easily calculated within the stan¬ 
dard WKB approximation. The total cross section follows from multiplying these 
eigenchannel cross sections by the weight factors indicated in Eq. (25). The calcula¬ 
tion of the weight factors (n\U\£i) within the IBM is straightforward and is given in 
Refs. |1|, g, to which the reader is referred for further details. 

Comparison of the linear and all orders coupling calculations shows that both 
can describe fusion cross sections equally well. However there are important differ¬ 
ences between the two approaches when one considers barrier distributions, and only 
through the latter, can one obtain a good representation of the data [12 . Another 


advantage of the all orders coupling is that the strengths used are very close to the 
deformation parameters and it is easier to establish a contact with the geometri¬ 
cal model. (In contrast, rather large strengths are required in the linear coupling 
approximation.) 









3 A Systematic Study of Subbarrier Fusion in Rare- 
Earth Nuclei 

We carried out a systematic study of subbarrier fusion data, accumulated over 
the last few years, for rare-earth nuclei JTB|. To illustrate the quality of this global 
fit, in Fig. 1 we compare our theoretical result for the 16 0 + 148 Sm, 16 0 + 154 Sm, 
and 16 0 + 186 W systems with the data from the Australian National University. 
There is excellent agreement between the IBM-based model and the data, as the 
target nucleus changes from vibrational ( 148 Sm) to deformed with positive ( 154 Sm) 
and negative ( 186 W) hexadecapole moments. 

Having extracted a consistent set of parameters, we can predict fusion cross sec¬ 
tions and barrier distributions in other rare-earth nuclei. In this respect, the transi¬ 
tional Os-Pt region is of particular interest. One expects that barrier distributions 
for transitional nuclei exhibit sharp changes due to shape-phase transitions, contrary 
to vibrational and rotational nuclei where the cross sections increase smoothly with 
increasing deformation or mass number [H|. An accurate measurement of subbarrier 
fusion cross sections for transitional nuclei such as Pt and Os provides a sensitive 
test for competing models in this region and possibly point to new directions in re¬ 
search. A subsequent measurement of the fusion cross sections for 40 C++ 194 Pt and 
40 C++ 194 Os systems confirmed our predictions [|21 . 


It is also possible to generalize the previous formalism to include arbitrary kinds 
of bosons in the target nucleus and investigate whether g bosons have any discernible 
effects on subbarrier fusion reactions. We found that |U| except for slight differ¬ 


ences in the barrier distributions (which can be made even smaller by fine tuning the 
coupling strengths), there are no visible differences between the sd and sdg model 
predictions. The similarity of the results implies that subbarrier fusion probes the 
overall coupling strength in nuclei, but otherwise is not sensitive to the details of the 
nuclear wavefuuctions. In this sense, subbarrier fusion is in the same category as 
other static quantities (energy levels, electromagnetic transition rates), and does not 
seem to constitute a dynamic probe of nuclei, in contrast to proton scattering 

Another experimental test of our model would be to study angular momentum dis¬ 
tributions in subbarrier fusion, which can be determined reasonably accurately from 
the gamma-ray multiplicities data or from the measurement of relative populations 
of the ground and isomeric states in evaporation residues 0. Vandenbosch and his 
collaborators at the University of Washington, Seattle devoted considerable time to 
the measurement of angular momentum distributions j6|. 191. In Fig. 2, we present 


angular momentum distributions calculated using our model for a number of systems 
0. We should emphasize that these angular momentum distributions are not di¬ 


rect experimental observables. However there are relations of rather general nature 


relating angular momentum distributions to the moments of fusion cross sections |22 


which can be used to assess the consistency of the direct measurements of the fusion 
cross sections with the angular momentum distributions extracted from gamma-ray 
multiplicities. 
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Figure Captions 

Fig. 1. Comparison of fusion cross sections and barrier distributions predicted 
using the Interacting Boson Model with the data as described in the text. The 
parameters are given in Ref. |K| . 


Fig. 2. Comparison of angular momentum distributions predicted using the In¬ 
teracting Boson Model with the data. The parameters are given in Ref. 
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